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We consider super-gravity models in which the lightest supersymmetric particle (LSP) is a stable 
gravitino. The next-to-lightest supersymmetric particle (NLSP) freezes out with its thermal relic 
density and then decays after (10 5 — 10 10 ) sec, injecting high-energy photons into the cosmic plasma. 
These photons heat up the electron plasma which then thermalizes with the cosmic microwave back- 
ground (CMB) via Compton scattering, bremsstrahlung and double-Compton scattering. Contrary 
to previous studies which assume instantaneous energy injection, we solve the full kinetic equation 
for the photon number density with a source term describing the decay of the NLSP. This source 
term is based on the requirement that the injected energy be almost instantaneously redistributed by 
Compton scattering, hence leading to a time-dependent chemical potential. We investigate the case 
of a stau NLSP and determine the constraints on the gravitino and stau masses from observations 
of the CMB spectrum by assuming that all gravitino LSPs come from stau NLSP decays. Unlike 
the analytical approximations, we find that there may be a stau mass below which the constraint 
from the CMB spectrum vanishes. 

PACS numbers: 95.35.+d, 98.80.Es, 98.80.Cq 



I. INTRODUCTION 

Supersymmetry provides mainly two compelling can- 
didates for cold dark matter: either the gravitino or the 
neutralino, depending on which one is the LSP. If we re- 
quire that R-parity be conserved, the NLSP decays into 
the stable LSP and releases energy in standard model 
particles. At leading order, these late decays are two- 
body and the accompanying energy is mainly electromag- 
netic. 

Part of the electromagnetic release is transferred to the 
cosmic microwave background radiation. The CMB then 
re-thermalizes through three relevant processes: Comp- 
ton scattering (7+e — > 7+e), double-Compton scattering 
(7 + e — * 7+7+ e) and bremsstrahlung (e+p — + e+p+7). 
The energy injection may distort the CMB, depending on 
the redshift at which it occurs and on the various time 
scales of the processes. Early NLSP decays can be fully 
thermalized, whereas distortions caused by injection from 
late decays cannot. Varying the NLSP and LSP masses, 
it is possible to control both the NLSP lifetime and the 
energy injected in the CMB. 

The CMB is not only very isotropic, but it also has 
a very precise Planck spectrum. The FIRAS instrument 
aboard the COBE satellite constrains the deviation from 
a perfect blackbody spectrum in terms of a few numbers. 
Important for this work is the limit for the chemical po- 
tential 



\H\ < 9 X 1(T 5 . 

This bound has been used to derive limits for the en- 
ergy released by NLSP decays as well as the NLSP life- 



time. Observational constraints on the CMB spectrum 
can be translated into bounds for the stau NLSP and 
gravitino LSP masses. It should be mentioned that we 
explicitly assume that all gravitinos present in the uni- 
verse are produced by stau decays. Some models suggest 
that gravitino LSPs may also be produced by scatter- 
ing interactions after reheating, leading to less stringent 
constraints. Recently, several papers have employed an 
analytic approximation to determine these limits 0-0-[l 
The approximation used in 121- [ii lif is derived in Ref. |5j . 
This analytical result turns out to provide the most strin- 
gent limit on the gravitino dark matter model in some 
range for the NLSP and LSP masses. This prompted us 
to repeat the calculation numerically. 

We find that the bounds for a chemical potential of 
fi < 9 x 10 given by |5J is a good approximation only for 
stau masses above 500 GeV. Below this mass, our bounds 
are less stringent and even disappear for staus lighter 
than 100 GeV. We also consider the limit /i < 10~ 5 and 
find that our results suggest lighter gravitinos or equiv- 
alently shorter stau lifetimes. Finally, we consider an 
upper bound on the chemical potential of 2 x 10 -6 as 
planned to be achieved in the DIMES experiment [||. 
We find that, if DIMES does not see a chemical potential, 
fi < 2 x 10~ 6 , gravitinos cannot significantly contribute 
to the dark matter if supersymmetry breaking is gravity 
mediated. 

After recalling some properties of the NLSP in the next 
section, we write down the full kinetic equation for the 
evolution of the photon one-particle distribution function 
in Section ITTT1 We express it as an evolution equation for 
a frequency-dependent chemical potential. This allows 
us to determine the chemical potential for a given point 



2 



in the (tonlsp, ™-lsp) and (m N LSP, tnlsp) plane of the 
NLSP and to draw exclusion plots in Section II VI We 
summarize our conclusions in Section Ivl 



II. NLSP PROPERTIES 

We investigate a super-gravity model with a gravitino 
LSP and a stau NLSP. We assume that the NLSP freezes 
out with thermal relic density and decays after a time 
determined by both its mass and the LSP mass. Due 
to the suppression of non-photonic decay channels, the 
branching ratio for decays to photons is set to be equal 
to one. 

Gravitino LSPs are produced through NLSP decays 
NLSP -► LSP + SMP, where SMP are standard model 
particles. Using the standard N = 1 super-gravity La- 
grangian, the rates for the various decay channels of the 
NLSP can be calculated. 



A. Slepton NLSP 



As long as the staus do not decay, their time-dependent 
number density can be expressed as 



n f (i) 0.26 m 



-3 



GcV 



T 



1 3 



Kelvin 



(3) 



The fact that the final gravitino density, fl^h 2 = 
(mQ/m,f)VLfh 2 , is bounded by observations, Vigh 2 < 0.14 
yj together with Eq. (J2J) implies an upper bound for the 
gravitino mass as a function of the stau mass: 



"!<L< .7xl0 6 ^. 



GeV 



(4) 



III. 



EVOLUTION EQUATION FOR THE 
PHOTON NUMBER DENSITY 



The decay of unstable particles into photons during 
the early stages of the universe can lead to distortions 
in the CMB. Depending on the redshift at which energy 
is injected, this may leave a measurable imprint of the 
early decays. This is the process which we now analyze 
in detail. 



We assume a gravity-mediated supersymmetry break- 
ing model, where the gravitino LSP mass is of the or- 
der 10 2 - 10 4 GeV and the stau NLSP lifetime of the 
order 10 4 - 10 10 sec. In gauge-mediated supersymmetry 
breaking models the gravitino is also the LSP but is much 
lighter (m^ < keV), resulting in a much shorter NLSP 
lifetime. In such models, the present CMB constraints 
do therefore not apply. 

The width for the decay of any sfermion / to a gravi- 
tino G for a negligible fermion mass is given by 



r(/->/<5) 



i 



48ttA/ 2 m 2 

Cjj 



(1) 



where M* = (8ttGn )~ 1 / 2 is the reduced Planck mass. 

Stau NLSPs decay to taus and gravitinos, which then 
decay to e, p,, ir°, and v. As mentioned in Ref. 
the electromagnetic energy produced in r decays varies 
between e^M ^E T and ej^f — E T . If not specified, a 
value of £em = 0.8_E r will be assumed throughout of this 
paper. 



B. Thermal relic density 

Using the thermal relic density of the right-handed 
slepton NLSPs determined in Ref.[| and the thermally- 
averaged cross section from Ref. y|, the stau relic abun- 
dance is given by 



ni h h 2 



0.2 



TeV 



(2) 



A. Energy injection by particle decays 

Energy injection resulting from NLSP decays heats the 
electrons, leading to a ratio T e /Tcmb larger than one, 
where T e is the electron temperature and Tomb is the 
temperature of the CMB. Due to the tight coupling be- 
tween electrons and photons during the early stages of 
the universe, the energy surplus of the electrons is redis- 
tributed among the photons, distorting the CMB pho- 
ton distribution from a blackbody spectrum. Assuming 
that the energy transfer between electrons and photons 
results in a Bose-Einstein spectrum with a frequency- 
independent chemical potential, it is possible to relate 
this resulting chemical potential to the number and en- 
ergy density of the injected photons and electrons. 

Following the analysis done in we can write the 
energy in a Bose-Einstein distribution as 



Pbe 



4(7 SB T 4 [ 1 



90C(3) 

_4 



(5) 



where we assume a small chemical potential. The number 
density is given by 



™BE 



2C(3)^ 3 



Ti 1 



6C(3) 



(0) 



Here £ denotes the Riemann (- function (see Ref. jlO|). 
Furthermore, due to energy conservation, we know that 
the energy density may also be written as 



PBE — PP + Pdocay, 



(7) 



where pdocay is the energy density injected by the NLSP 
decays and pp is the density of the CMB photons. This 
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equation is only valid if the injected photons are re- 
distributed in a negligible amount of time compared to 
the time scales of double- Compton and bremsstrahlung. 
However, this does not hold for the low-frequency spec- 
trum, where the photon-creating processes dominate. 

More precisely, the injected energy density is given by 
the following differential equation 



dp, 



'decay 



'NLSP 



'LSP 



dt 



2TONLSP 



-t/r 



.a 

^ Pdecay > 



(8) 

where r is the lifetime of the NLSP. Due to the fact that 
tau decays also produce several neutrinos, the right hand 
side of Eq. JSJ) has been multiplied by a factor e describing 
the ratio of the injected energy to the total energy. As 
pointed out in Sec. Ill Al e may have a value between 0.3 
and 1; we set e = 0.8. In order to solve Eq. we 
integrate both sides from t m to t and obtain 



Pdccay 



"NLSP 




where Erf is the error function as defined in Simi- 
larly, the photon number density is given by 



7J B E = ™P + ride 



(10) 



where "decay is the injected photon number density given 
by the following differential equation 



dn, 



decay 



= A^nNLSp(*)- 



-t/T 



1 — 

<J "decay- 



(11) 
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FIG. 1: Solution of the system of equations 1131 . Top: Time 
evolution of the chemical potential for fixed Ny at m? = 200 
GeV and = 50 GeV, corresponding to tnlsp ~ 6 x 10 6 sec 
(blue solid line) and vtif — 300 GeV and = 50 GeV, 
corresponding to tnlsp ~ 7 x 10 5 sec (green dashed line). A 
value of iV 7 = 10 4 has been assumed. Bottom: The chemical 
potential at a time t — 10 8 sec is shown as a function of iV 7 . 
Both lines are given by the same masses as above. 



Here Ny is the number of photons per stau decay injected 
into the spectrum. The solution of Eq. I|ll|) is given by 



"decay = NyUNLSP (*in) (~) 



(12) 



Inserting Eq. and Eq. Q into Eq. J5J , as well insert- 
ing l(T2|> and Eq. (fTUf> into Eq. © we find the relations: 



"decay 

n p 

Pdccay 
PP 



T, 



6((3) 
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(13) 



Up to now, we have assumed that the energy injection 
caused by stau decays was instantaneously converted in 
a chemical potential through Compton scattering. How- 
ever, there are also two other processes that do not con- 
serve the photon numbers: bremsstrahlung and double- 
Compton scattering. The influence of both processes is 
discussed in the next section. 

We shall use the value /ij n j obtained in this section as 
initial condition for the numerical solution of the Boltz- 
mann equation discussed below. 



B. Photon- matter interaction 



These equations cannot be solved simultaneously since 
there are three unknowns, /^inj, T e and Ny. However, it 
turns out that the chemical potential is independent of 
Ny up to an unreasonable photon number injection of 
w 10 7 photons per NLSP decay, as shown in the bottom 
panel of Fig. 2] The top panel shows the time evolution 
of the chemical potential for two different stau lifetimes. 



When the universe is more than a few minutes old, the 
coupling of CMB photons and matter is basically due 
to three processes: Compton scattering, double Comp- 
ton scattering and bremsstrahlung. If these processes 
are no longer very efficient, the spectrum can be dis- 
torted. Especially, if double Compton scattering and 
bremsstrahlung become weak, the photon number can 
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no longer be changed and energy injection into the CMB 
leads to a chemical potential. With this in mind, we 
parameterize a general distorted spectrum by a Bosc- 
Einstein distribution with a frequency-dependent dimen- 
sionless chemical potential /j,(x,t), 



n(x, t) 



exp(x + fj,(x, t)) — 1 



(14) 



where x = hv/T e is the dimensionless photon frequency. 
The collision terms in the Boltzmann equation for the 
three relevant processes (Compton scattering, double- 
Compton scattering and bremsstrahhmg) are studied in 
[TT| . The Boltzmann equation is given by 



dn 
~dt 



1 T e 1 d (dn 

yg /Tig X Q 

Qg(x) 1 



t 7e m e x 2 dx V dx ^ 



4- fj X ;-f» 3 

L/y g O Jj 



[l-n(e*-l)) 



1 4a / T P V 1 



^ ^ ^ t 7 e 37T \m e ) X 3 

x [l-n{e x -1)} J dxx 4 (l + n) 

(gZ+Wnj _ I) 2 dt 



(15) 



The first term describes Compton scattering, the sec- 
ond bremsstrahlung, the third term double Compton 
scattering and the fourth is the injection term given 
by the solution of Eq. Ijl3|l . We have introduced the 
Heaviside function 9 in the double-Compton scattering 
term, to take into account that it is active only for 
x < 1. The constant t ie — {n e OT)~ x is the Thom- 
son scattering time, Q — 2^/2ii(m e /T e ) 1 / 2 anBT~ 3 ~ 
1.7 x 10 _10 (MeV/T) 1 / 2 (T/T e ) 7 / 2 f7 B /i 2 and g[x) is the 
Gaunt factor. More details on the collision terms can be 
found in Refs. H HH or 

Since we expect a small value of the chemical potential, 
we can expand this equation to first order in /j. 



n(x,t) fa n (x,t) +/^(x,t)^-(a;,i) 
= ~z 7 _ l ) 



(e x - l) 5 



(16) 



The zeroth order is the equilibrium distribution, and 
the first order in \i describes the spectral distortion. 
The kinetic equations for the three relevant processes 
(Compton scattering, double Compton scattering and 
bremsstrahlung) [Tl| , then becomes a linear equation for 



the evolution of the chemical potential (// = d[i/dx), 







(e x - l) 2 dt 



/j,(x,t) 



2 T„ 



xe 



2x 



t je m e (e x - I) 4 



[(4 — 4cosha; + x sinh x)n'(x, t) — x(coshx — l)fi"(x,t)] 



H{x,t) 



+ [1 - 6(x - 1)] 



1 16n 3 a 



45 



-u(x,t) 



x 3 (e x — 1) 

dt 



(e x - 1)' 



(17) 



We have solved both systems of equations numerically 
and find consistent results. When interested in values 
/i ~ 10 -5 , we start at t m = 10 5 sec, but when we want 
to detect chemical potentials on the level of \i ~ 10~ 6 we 
have to start at t; n = 10 4 sec. 



C. Time evolution of the frequency-dependent 
chemical potential 

As shown in Ref. energy injected at a redshift 

higher than z ~ 10 7 , corresponding to a time t < 10 5 sec 
is fully thermalized. Furthermore, at decoupling time 
idec ~ 10 13 sec, the CMB spectrum is frozen in and does 
not evolve anymore apart from redshifting the photon 
momenta. However, as shown in Fig. [21 the photon- 
creating processes are unable to reduce a chemical po- 
tential already as early as t > 10 s sec, much before re- 
combination. 

We see from the top panel of Fig. [21 that, compared 
to the chemical potential of a Bose-Einstein spectrum 
(dashed blue line), double-Compton and bremsstrahlung 
significantly reduce the magnitude of the distortions from 
a blackbody spectrum: the chemical potential at late 
times has been reduced from 7.4 x 1CP 5 to 4.3 x 1CP 5 . The 
bottom panel of Fig.[2]shows the frequency dependence of 
fi evaluated at different times. The high-frequency range 
is dominated by Compton scattering. The chemical po- 
tential is constant above x > 1, describing a true Bose- 
Einstein spectrum. The low-energy spectrum is domi- 
nated by the photon-creating processes which can destroy 
the chemical potential and lead to a Planck spectrum be- 
low i<4x 10~ 3 at recombination time. 

It is clear that the later the energy injection, or equiv- 
alently the later the staus decay, the weaker are the 
photon-creating processes which would reduce the distor- 
tions. However, it should be kept in mind that our equa- 
tions are valid only if Compton scattering can achieve 
a Bose-Einstein spectrum. By requiring that Compton 
scattering be well active during stau decays, we can put 
an upper limit on the stau lifetime. Following the anal- 
ysis of |T3 |. a given spectrum can only relax to a Bose- 
Einstein spectrum before teE — 10 9 sec. Therefore, the 
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X 



FIG. 2: Solution of the full kinetic equation 115H for m? = 200 
GeV and uiq = 50 GeV, corresponding to tnlsp ~ 6x 10 6 sec. 
Top: The dashed (blue) line represents the chemical potential 
given by Eq. 1131 . the solid (red) line is the solution of Eq. (1151 
evaluated at x = 3. Bottom: Chemical potential as a function 
of x given by Eq. 1151 evaluated at t = 5 x 10 6 sec (dotted 
blue line), t = 10 7 sec (dashed green line) and t = idee (solid 
red line). 



accuracy of a solution of Eq. I|15|l for a stau lifetime longer 
than tnlsp ~ 5 x 10 s sec becomes questionable and un- 
trustworthy for tnlsp <L 10 9 sec. 

On the other hand, after freeze-out of Compton scat- 
tering, the injected energy resulting from stau decays 
cannot be scattered downward in frequency. But there 
are several other processes that could leave an imprint 
on the measurable CMB spectrum. For example, pho- 
tons produced during stau decays have an energy much 
greater than the electron mass and can create electron- 
positron pairs through the process 7 + 7 — * e + + e~. 
The rate of this process has a typical value of Tdp ~ 10 3 
sec -1 , provided E 1 > m 2 e /22T 13j. Compared to the 
Hubble rate H = l/2t, we see that this process plays 
an important role in heating up the electrons. However, 
the raise of the electron temperature does no longer af- 
fect the CMB spectrum when Compton scattering has 
already frozen out. 

We have also analyzed the fact that the true electron 
temperature T ei t ru e is not the same as the electron tem- 
perature T e obtained by solving Eq. (|13fl due to the influ- 



ence of bremsstrahlung and double-Compton scattering 
which reduce the chemical potential from the value jUi n j- 
Given the chemical potential fi at recombination time, 
we can calculate T e by inserting \x into one of the equa- 
tions of the system Ill3fl . We have found that T e true only 
differs from T e by ~ 10~ 4 , and that changing the electron 
temperature in Eq. I|15f) by such small amounts has no 
effect. (This is not so surprising, as the effect is of second 
order.) 




FIG. 3: Excluded and allowed regions in the (ttinlsp, tilsp) 
parameter space. The light shaded region in the upper part 
of the graphs, labeled 'G not LSP', (yellow) is not considered 
since in this part of parameter space jtinlsp < ttilsp- The 
shaded region in the right upper corner, labeled 'fig', (green) 
is excluded by the over-closure constraint fi ( =i/i 2 ;$0.14. The 
dark shaded region on the left, labeled 'LEP data' (red) is 
forbidden by LEP measurements [Tflj . The dashed line cor- 
responds to QqIi 2 = 0.12, the best fit value for the cold 
dark matter density from the WMAP experiment . The 
dark shaded region in the middle (magenta) is forbidden for a 
chemical potential limit of fi < 9 x 10 -5 (top) and fi < 10~ 5 
(bottom). The approximation of Ref. |5J excludes the entire 
region above the solid line. 
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FIG. 4: Excluded and allowed regions in the (ttinlsp, t) plane. 
The dashed line corresponds to Q^h 2 — 0.12, the best fit 
value for the cold dark matter density from the WMAP ex- 
periment [jj. The regions are labeled as in Fig. El In the top 
panel the limit on the chemical potential is /i < 9 x 10~ 5 while 
we require fj, < 10 -5 in the bottom panel. The approximation 
of Ref. |3 excludes the entire region above the solid line. 



IV. CMB CONSTRAINTS ON THE STAU AND 
GRAVITINO MASSES 

The FIRAS instrument aboard the COBE satellite has 
measured a temperature To = 2.725 ± 0.001 Kelvin 
and it was able to give an upper bound for the chem- 
ical potential H < 9 x 10~ 5 . This bound 
comes from measurements in the frequency range from 
2 to 600 GHz, corresponding to x = hv/T € [0.03, 10]. 
There are also some measurements at lower frequen- 
cies, but their accuracy is worse, leading to a lower 
bound on [i which is by at least an order of magnitude 
higher. To obtain good accuracy in the measured inter- 
val, we numerically compute the chemical potential for 
x E [10~ 4 ,15]. We require that the chemical potential 
be never higher than 9 x 10~ 5 within the experimental 
range x £ [0.03, 10]. Outside that range, [i may be larger 
(experiments do not rule out deviations outside the fre- 
quency range [0.5 GHz, 600 GHz]). 

A point in the (tonlsPj WLSp)-P ian e is considered to 




FIG. 5: Excluded and allowed regions in the (tunlsp, wlsp) 
plane for different values of the parameter e and a chemical 
potential (i = 9x 10~ 5 . The labeled regions are as in Fig. [3] 
The lower solid (blue) line represents the case e = 1, the 
dashed (red) line e = 0.8 and the dotted (black) line e = 0.3. 
The upper solid line denotes the limit on the lifetime, r > 
10 9 sec above which the CMB spectrum is not modified. 



satisfy the CMB observational bound if the magnitude of 
the chemical potential never trespasses the limit 9 x 10~ 5 
within and only within the frequency range [0.03, 10]. 
Due to the limitations explained in the previous section, 
not every point in the (tonlsp, wilsp) -plane can be cal- 
culated, but we expect the chemical potential to be much 
smaller than the experimental limit for points where our 
calculation cannot be trusted. 

An estimate of the chemical potential caused by an in- 
stantaneous energy injection is given in |5j. It is shown as 
a solid line in Figs.|31to[S] This approximation is not very 
precise for small distortions. While it is in good agree- 
ment with our results for staus heavier than 500 GeV, 
it does not take account of the fact that light staus do 
not inject enough energy to significantly distort the spec- 
trum. As shown in the upper panel of Fig. |3J the bound 
on the gravitino mass disappears for staus lighter than 
« 100 GeV. Moreover, due to the freeze-out of Compton 
scattering, we have introduced a limit on the gravitino 
masses corresponding to a stau lifetime 10 9 sec. All grav- 
itino masses leading to longer stau lifetimes are allowed. 

Our results match quite well with the approximation 
given in Q for fi < 10~ 5 , as shown in the bottom panel 
of Fig. |3 We find the same limit for light staus and a 
somewhat more stringent limit for heavy staus. However, 
when compared to the stau lifetime (see the bottom panel 
of Fig^J , our numerical results give a limit on the lifetime 
that is up to five times shorter than the one obtained by 
using the approximation. 

As mentioned in Sec. lIII A1 the injected energy depends 
on the energy going into neutrinos. This is described by 
the parameter e. Contrary to what is claimed in 
our results depend significantly on this parameter. We 
considered the two extreme cases e = 0.3 and e = 1 
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and the usual case e — 0.8 (Fig. |SJ). The case e = 0.3 
is much less constraining than e = 0.8, even completely 
disappearing for staus lighter than 140 GeV. On the other 
hand, the cases e = 1 and e = 0.8 match well down to 
rrif ~ 200GeV below which the former becomes more 
stringent. 

Future missions like the Absolute Radiometer 
for Cosmology, Astrophysics Diffuse Emission (AR- 
CADE) 0, Il7l or the Diffuse Microwave Emission Sur- 
vey (DIMES) [_J experiments may improve sensitivities 
in the poorly studied centimeter-wavelength band, im- 
proving the limit on the chemical potential to about 
< 2 x 10~ 6 . In our model, if neither DIMES nor AR- 
CADE is able to measure distortions of the CMB, grav- 
itinos could only contribute to the missing dark matter if 
t nlsp > 10 9 sec or tnlsp £ 2x 10 5 sec (see Fig.®. How- 
ever, combining our results with other constraints |18| . 
we find that gravitinos could not significantly contribute 
to the dark matter for such a bound. 




FIG. 6: Excluded and allowed regions in the (mNLSP, wlsp) 
plane (top panel) and in the (ttznlsp , t) plane (bottom panel) 
for a limit on the chemical potential of fj, < 2 x 10~ 6 . The 
labeled regions are like in Fig. 01 



V. SUMMARY 



We have studied the effect on the CMB from stau 
NLSP decays into gravitino LSPs, assuming that the 
staus freeze out with their thermal relic density. We 
have numerically solved the kinetic equation for the pho- 
ton number density with non-instantaneous energy in- 
jection. We have found that our numerical results are in 
good accordance with the analytical approximation [_j 
for the induced chemical potential fi < 9 x 10 -5 if the 
stau is heavier than 500GeV, but differs considerably for 
lighter stau masses. For light staus the constraints are 
weaker and even disappear for m f < lOOGeV. On the 
other hand, the approximation underestimates the limits 
for stronger constraints given by fi < 10~ 5 or even more 
for /i < 2 x 10~ 6 . This limit, which could be achieved 
in planned experiments [_i llfil| , together with other con- 
traints would completely exclude the gravitino as 
dark matter candidate in models with gravity-mediated 
supersymmctry breaking. However, allowing a gravitino 
production after reheating leads to less stringent con- 
straints than our results. We also found that the results 
depend sensitively on the energy injection parameter e. 
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